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We examine the outer Galactic HI disk for deviations from the 6 = 0° plane by 
constructing maps of disk surface density, mean height, and thickness. We find 
that the Galactic warp is well described by a vertical offset plus two Fourier 
modes of frequency 1 and 2, all of which grow with Galactocentric radius. A 
perturbation theory calculation demonstrates that the tidal influence of the 
Magellanic Clouds creates distortions in the dark matter halo, which combine 
to produce disk warp amplitudes comparable to the observations. We also 
construct a detailed map of the perturbed surface density of HI in the outer 
disk demonstrating that the Galaxy is a non-axisymmetric multi-armed spiral. 
Overdensities in the surface density are coincident with regions of reduced gas 
thickness. 



1 Method 

We use the 21cni Leiden/ Argentine/Bonn (LAB) data (Kalberla et al., 2005; 
Hartmann & Burton, 1997; Bajaja et al., 2005; Arnal et al., 2000) to conduct 
a quantitative study of the warp and spiral structure. We used the Banning 
smoothed data, which have a velocity resolution of 1.9 km s~^, and restricted 
our sample to |6| < 30°. We apply a median filter to the data to remove 
angularly small, bright features that are clearly not associated with the disk. 
To create a map of the Galaxy from the data, we must assume a rotation 
structure for the disk. Using a flat rotation curve with circular orbits results 
in a large asymmetry between the surface densities at Galactic longitudes 
on either side of £ = 0° and t = 180°. To eliminate this discontinuity, we 
assume that the gas is traveling on elliptical streamlines configured such that 
the surface densities are smooth across the i = 180° line (see Levine et al. 
(2006a)). We exclude aU points that lie within 345° < € < 15° or 165° < 
i < 195°. Points in these two wedges have velocities along the line of sight 
that are too small with respect to their random velocities to establish reliable 
distances. 
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With this rotation structure, we transform the data to a three dimensional 
density grid, p{R,(j),z), where R is Galactoccntric radius and z is height off 
the plane. The Galactocentric azimuth </> is defined such that it converges 
with Galactic longitude I at large R. We use a dispersion filter on p to remove 
extended clouds near the disk or spurs that have split off from the disk. 
Finally, we calculate the surface density S and mean height /i as a function of 
position [R, </>) in the disk. The half-thickness Th is calculated from p before 
the dispersion filter is applied, following the procedure described in Henderson 
et al. (1982). 

2 Warping 

2.1 Modes 

A Lomb periodogram analysis of each radius ring reveals that the power in 
each of the frequency m = 0, 1, and 2 modes is consistently larger than that 
in any other mode for i? > 20 kpc. Accordingly, we characterize the warp by 
an offset in the z direction, plus two Fourier modes with frequency m = 1 and 
2. We fit each ring with the function 

W{(j)) = Wq + Wi sin((^ - (^i) + W2 sm{2(l> - (l>2); (1) 

this is similar to a fit described in Binney & Merrificld (1998). Each of the 
three amplitudes Wj and two phases cjji in this fit is a function of radius, 
because we fit each radius ring independently. This three component fit does 
a good job of reproducing the large-scale features in the mean height map 
(Fig. 1). 

The three amplitude parameters each increase monotonically, with the 
m = mode possibly reaching an asymptotic value near the far end of our 
radius range. At i? « 11 kpc, the m = 1 mode dominates the shape of the 
warp; the other two modes do not become important until i? w 15 kpc. There 
is little evidence for precession in the lines of maxima for the rn = 1 and 2 
modes, and the line of maxima for the m = 1 mode is roughly aligned with 
one of the lines of maxima of the m = 2 mode. 

2.2 Mechanism 

Many efforts have been directed toward understanding the warp on a the- 
oretical basis. Bending modes have long been suspected as the mechanism 
creating and maintaining the warp (Hunter & Toomre, 1969; Sparke & Caser- 
tano, 1988: Sparke, 1995, and many others), but several other mechanisms 
have been suggested (Kuijken & Garcia-Ruiz, 2001). Gravitational interac- 
tion with satellites such as the Magellanic Clouds has been a long-standing 
possibility (Burke, 1957; Weinberg, 1998), but there is a debate over whether 
tidal effects are strong enough to produce the observed effect (Kerr, 1957). 
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Fig. 1. Left panel: h{R,4>)- Right panel: The fit to the warp is plotted, along with 
the lines of maximum amplitude for the m=l (dotted) and m=2 (dashed) modes. 
These lines are marked Wi and W2 respectively. The white contour line denotes a 
height of kpc; black lines mark the 2, 4, and 6 kpc elevations. The colorbar is 
marked in kpc. 



We perform a linear perturbation analysis to model the disk's response to 
the Clouds (Weinberg & Blitz, 2006). We allow the disk to feel the tidal 
field from the Clouds directly as well as the force from the dark-matter 
halo wake excited by the Clouds. We adopted the LMC mass from West- 
erlund (1997), the radial velocity and proper motion from Kallivayalil et al. 
(2006), and the distance modulus from Freedman ct al. (2001). The warp 
is a very dynamic structure based on the temporal evolution of the model. 
This can be seen in the AVI file of the simulations which can be found 
at http:/ /www.astro.umass.edu/^weinberg/lmc. Rather than a static 
structure that might be expected for a warp in response to a triaxial halo, a 
warp that results from the Clouds is continuously changing shape because of 
the varying amplitudes and phases of the various modes. As in the data, the 
m = 1 mode in the calculations is the strongest and increases nearly linearly 
out to the edge of the disk. The calculations also show a weak response of 
m = 0, 2 modes out to about 15 kpc, and then increasing nearly linearly, but 
with approximately the same amplitude for both. All of the major features of 
the warp are reproduced by the calculations, though there are some relatively 
minor differences which are likely to be due to uncertainties in the modeling. 
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3 Spiral Structure 

Mapping the Milky Way's spiral structure is traditionally difficult because 
the Sun is imbedded in the Galactic disk; absorption by dust renders optical 
methods ineffective at distances larger than a few kpc. Radio lines like the 21 
cm transition are not affected by this absorption, and are therefore well-suited 
to looking through the disk. In §1, we described the calculation of E{R,<j)); 
we now perform a modified unsharp masking procedure on this map (Levine 
et al., 2006b). By subtracting a blurred copy from the original image, this 
technique emphasizes low contrast features in both bright and dim regions. In 
previous maps (Henderson ct al., 1982) spiral arms were difficult to pick out 
on top of the global falloff of surface density with radius. 

We construct a ratio n{R, cj)) between the surface density and the local 
median surface density. This is a dimcnsionlcss quantity, and therefore is a 
direct measure of the strength of the surface density perturbations (Fig. 2). 
The vast majority of points have values in the range 0.6 to 1.8, implying that 
the arm-to- intcrarm surface density ratio is about 3. Several long spiral arms 
appear clearly on the map, but the overall structure does not have the origin 
reflection symmetry of a "grand-design" spiral. Easily identified spiral arms 
cover a larger area in the south than in the north; there arc three arms in the 
southern half of the diagram. Each of these arms has already been detected 
in previous maps out to 17 kpc for the two arms closer to the anticenter 
(Henderson et al., 1982) and out to 24 kpc for the 'outer' arm (McClure- 
Griffiths et al., 2004). 

We also perform the modified unsharp masking on the thickness Th{R, (p); 
spiral structure is evident in this map as well. Wc plot the 77 = 1.1 ovcrdcnsity 
contour from the surface density perturbation map on top of the thickness 
perturbation map (Fig. 3), showing that there is a good match between the 
arm positions as calculated from the surface density and the thickness; the 
thickness of the HI layer is smaller in the arms than in the rest of the disk, as 
suggested previously (Henderson et al., 1982). 
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Fig. 2. A contour plot of the ratio of the surface density to the local median surface 
density n{R,4>). Shaded regions are overdense compared to the local median. The 
thick solid contour marks the line 11 = 1.1, while the thinner contour marks the line 
n = 0.7. The values of U for the different contour levels are given by the colorbar. 
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Fig. 3. A contour plot of the perturbations in the gas thickness, r{R,4>). Shaded 
regions have reduced thickness compared to the local median. The perturbation 
levels for the different contour levels are given by the colorbar. The solid contour 
marks the line 77 — 1.1, the same as in Fig. 2, to show the alignment of the overdense 
surface densities with the thinner gas regions. 
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